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1. Introduction 



The supersymmetric M = 2 theories in two and four dimensions are quite interesting 
both for physicists and mathematicians. Their low-energy dynamics is rich, yet the super- 
symmetry is large enough to impose rather stringent constraints on the effective action. 

Prepotential of effective theory 

In quantum field theory one is interested in the low-energy effective action. The extended 
supersymmetry implies that all the terms with at most two derivatives and at most four 
fermions in the effective action are related, and can be expressed through derivatives of 
a single holomorphic quantity, a locally defined function jF(a) on the moduli space M. of 
vacua, the so-called prepotential. 

In the recent years the exact calculations of the prepotential were performed, which 
utilized the fact that only the one-loop perturbative corrections about (an arbitrary) instan- 
ton solution summed over the instanton sectors, contribute to the prepotential. Actually, 
the prepotential itself is not so easy to calculate. However, it turns out that a certain 
"quantum corrected" prepotential 



has a simple expression in terms of the gauge theory instantons. The leading term J- (a) = 
Fo(a) is the prepotential of the low-energy effective theory. The higher order terms, F g (a), 
describe the coupling of the theory to the Af = 2 supergravity multiplet. In the string 



string amplitudes. 

The prepotential F and the higher genus corrections F g are special in the sense that 
they determine the terms in the Lagrangian of effective theory, given by the integrals over 
half of the superspace, the so-called generalised F-terms. Formally one can consider the 
similar F-terms in the ultraviolet theory. The original microscopic Lagrangian 



oc 




(1.1) 



9=0 



theory realisations of the Af = 2 theory the terms F g (a) are computed by the genus g 




(1.2) 



can be deformed by adding the operators 




(1.3) 
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where 

<f> = <$> + 6ip + 96F- + ... (1.4) 
is the vector superfield, and it is often convenient to work with the generating function 

fc>0 

The bare couplings to, To are adjusted so as to produce the finite couplings in the effective 
Lagrangian, obtained by integrating out the high frequency modes. 

Renormalisation group flow 

We can reformulate our problem above as the calculation of the effect of the Wilsonian 
renormalisation group flow on the prepotential of the theory. One starts with the theory 
with Af = 2 supersymmetry, which is determined by the ultraviolet prepotential J-'vv, 
perturbed by arbitrary powers of the holomorphic operators 

^uv = 5(70 + *i) tr $ 2 + ^t fc - ' ® + (1.6) 



fc + 1 

fc>0 



and quadratic J-'jjy, i.e. 



— 1 —2 
•^uv = oT tr $ 



Then one integrates out the fast modes, i.e. the perturbative fluctuations with momenta 
above certain scale 11 as well as the non-perturbative modes, e.g. instantons (and fluc- 
tuations around them) of all sizes smaller then yT 1 . The resulting effective theory has a 
derivative expansion in the powers j^. The leading terms in the expansion are all deter- 
mined, thanks to the Af = 2 supersymmetry, by the effective prepotential J r ( / u). As \i is 
lowered all the way down to zero, we arrive at the infrared prepotenial 

^UV ► -^IR 

The supersymmetry considerations suggest that the renormalisation flows of T and T 
proceed more or less independently from each other. Thus one can simplify the problem 
by taking the limit, To — > «oo, while J^uv kept fixed. In this limit the path integral is 
dominated by the gauge instantons. The setup of [2] allows to evaluate their contribution, 
as well as the contribution of the fluctuations around the instantons, exactly. The price 
one pays is the introduction of extra parameters into the problem, some sort of the infrared 
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cutoff, which we denote by h~ 2 , since it appears to be a parameter of the loop expansion 
in dual topological string theory [1]. As we send % — > 0, the infrared cutoff is removed, 
and the prepotential is recovered as the extensive part of the free energy, cf. (1.1). For 
details the reader is invited to consult [2] [1] [3] . 

Prepotential and symplectic geometry 

The goal of this paper is to extract the generalised prepotential, JF IR = J 7 (a, t) as a function 
of the moduli of the vacua a, and the higher Casimir couplings t. It has been proposed 
in [4] (see also [5] and references therein) that the answer is given by tau-function of 
a quasiclassical or universal Whitham integrable hierarchy [6] and below we derive this 
hierarchy directly from the results of instanton calculus. 

For fixed t the prepotential defines a generating function of a Lagrangian submanifold 
C t in the complex symplectic vector space C^, invariant under the action of a certain 
discrete subgroup V of Sp(2iV, Z), the group of electric-magnetic dualities [7]. The defor- 
mations of such submanifolds are constrained by the considerations of duality, however at 
each order of deformation beyond the first one there are ambiguities: 



^(a, t) = F(a, 0) + ^ tkUk+i 

k>0 



(1.7) 



where Uk ~ (tr$ fe ) are T-invariant functions on £ , i.e. the polynomials of the coefficients 
of Seiberg-Witten curve, Cu (it) are some polynomials of it, called Losev-Shatashvili polyno- 
mials, and the precise definition of the theta constant and the theta characteristics in (1.7) 
is at the moment immaterial. While the choice of Uk is more or less a matter of convention, 
the choice of Losev-Shatashvili polynomials Cki(u) is physically important. The appear- 
ance of theta functions in (1.7) is dictated by the T-invariance and the extra terms Cki(u) 
are constrained by the degree considerations [7], which predict that for small k,l, they 
vanish. However, to determine them one needs a microscopic theory, and below we show 
how a particular choice of microscopic theory determines all contact terms Cki (u) in terms 
of quasiclassical integrable hierarchy. In particular, we shall see that Losev-Shatashvili 
polynomials Cki(u) indeed vanish for k,l < N, i.e. exactly for the Casimirs which deform 
the Seiberg-Witten geometry for the U(N) gauge theory. This distinguishes the quasi- 
classical hierarchy derived in this paper as directly coming from the microscopic instanton 
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theory, compared to the previous attempts, proposed in [8], based on the analogies with 
Landau- Ginzburg models. 



Ensembles of partitions 

In particular, for the so-called noncommutative U(l) theory (which supports instantons, 
[9]), or the theory on a single D3 brane in the background, which preserves only sixteen 
supercharges (so that the theory on the brane has only eight supercharges), the instanton 
partition function Z(a, h, t), t = (t\, t^-, ■ ■ ■), can be shown to be given by the sum over the 
Young diagrams, i.e. over the partitions [2][1][3]: 

7i ± %\ ST m A 1 ch fc+ i(a,A) 

Z(a,t,n) = > 7^-— exp > tk ; : (1.8) 

Y (~ h )' 1 h too k + 1 

(we shall remind the relevant notions of the theory of partitions in the Appendix). This 
theory can also be realised at a special point on the moduli space of U(N) gauge theory 
with 2N — 2 fundamental hypermultiplets' . If the theory has the gauge group U(N), e.g. 
it is realised on the stack of N fractional D3 branes, the corresponding partition function 
is given by the generalisation of (1.8): 

Z(3,t,ft) = ^ ert (a,t,n)^(m(a,A,n)) 2 (-l)l X l e xp ^ £ ^fc + 1 ^ (1-9) 

A k >° 

where m(a, A, %) is the U(N) generalisation of Plancherel measure [2] [3] and Z pert (a, t, %) 
is the perturbative partition function. 

Limit shape and Toda chain 

We shall evaluate the sums like (1.8) by the analogue of the saddle point method. The limit 
h — > of the sum (1.8) is dominated by a partition A*, of a large size ~ h~ 2 . The shape 
of the Young diagram of this partition, the so-called limit shape, is found by extremizing 
the effective energy functional. When all but the first coupling t\ are set to zero, this 
limit shape is the celebrated "arcsin law" curve, found by Vershik-Kerov and Logan-Shepp 



t In a certain sense, all the instanton contributions in the "U(l)" theory are the artefacts of 
the imbedding of the theory into the theory with rich ultraviolet structure. In the terminology of 
[10] these are "freckled instantons" and their contribution must be subtracted by the appropriate 
"mirror map". See also [11] for the related and more detailed discussion, and [12] for the related 
M = 1 considerations. 
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[13] [14] [15] [16]. Our main claim for the rank N = 1 is that the evolution of the limit 
shape under the higher Casimirs is governed by the quasiclassical Toda hierarchy, and we 
present the corresponding solution explicitly. This result, which we prove in the section 
3, gives a further confirmation of the claim of the previous work [1] of A. S. Losev and 
the authors, that the BPS sector of the rank one Af = 2 gauge theory is equivalent to the 
stationary sector of the CP 1 Gromov-Witten theory. The dispersionless Toda hierarchy is 
well-known to describe the genus zero part of that theory (see, e.g. [17]). 

In the section 4 we describe the non-abelian theory. We find that the corresponding 
limit shape is described by the Krichever quasiclassical tau-function [6] associated with a 
family of hyperelliptic curves with two marked points 

TV 

y 2 = ]J(z - x+(a, t))(z - xY(a, t)) (1.10) 
i=i 

Interestingly enough, for nonvanishing t this family extends outside the family of Seiberg- 
Witten curves. The answer is encoded in (1.10) and a particular (1, 0) -differential <i£, 

rfs = a(z;a,t)dz 

y 

holomorphic outside the points P± where z = oo, which obeys certain normalisation condi- 
tions. We demonstrate how the instanton corrections (for nonvanishing t) can be extracted 
from the generalised Seiberg-Witten geometry, and derive equation (1.7) from the "first 
principles" . In the section 5 discuss the relation of our approach with the Eguchi-Yang 
matrix integral and propose a possible way to extend this relation beyond the quasiclassical 
theory. 



6 



2. Gauge theory partition function 

We study M = 2 gauge theory in the self-dual Q-background. The setup and the rel- 
evant physics are reviewed in [2] [1] [3] so we just briefly discuss it here. The gauge theory 
path integral is saturated by instantons, and Af = 2 supersymmetry cancels the contribu- 
tion of fluctuations, so that effectively one has to integrate " unities" over the corresponding 
moduli spaces of the instanton solutions. These moduli spaces can be described by the 
ADHM construction, and the corresponding integrals can be computed via equivariant lo- 
calisation technique in nontrivial O-background [2]. As a result, the gauge theory partition 
function can be presented effectively in one of the following ways: 
1) The sum over instantons can be interpreted as a Van-der-Vaals gas in one dimensions: 

the integrals over the instanton moduli space of charge k can be reduced, via ADHM 

construction and equivariant localisation a la [18] to the grand canonical ensemble of 

a one-dimensional gas of particles 



Z(a,t,ei, 2 ) =ZP ert (a,t,e 1>2 ) 



x 



X V — ( 61+62 1 <T) I I do, exi) I L (O, ) I x 



j[ d0 je xp(--i-*7(0j)) 
JJexp (-J-Vfa-tj)) 



x 

KJ 

where U(x) = t(x) + t(x + e\ + e 2 ) — t(x + ei) — t(x + e 2 ) + 

+ eif2 ^ log ((x - aif - - (ei + e 2 ) 2 J 
l=i ^ ' 



(2.1) 



V(x) = -eie 2 log 



x 2 (x 2 - (ei + e 2 ) 2 ) 



(x 2 — el) (x 2 — e|) 

interacting via pair-wise Van-der-Vaals kind of potential, in the presence of (for 
U(N) theory) sources. Here ei )2 are parameters of the O-background, which can be, 
after calculating integrals in (2.1) via residues, put to be e\ = — e 2 = %. 
2) The gauge theory partition function can be shown to be equal to the partition function 
of the statistical model of random partitions. For U(N) gauge theory one deals with 
the ensemble of A-tuples of partitions A = (Ai, . . . , Ajv), and the partition function of 
the gauge theory is [2] [3] : 

Z(a, t, h) = (a, t, h) ■ m2 («> ^ h ) ("I) |X| exp^ £ t k chk +^ (2.2) 

x k=l 
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3) The gauge theory partition function can be also written as a path integral in the 
theory of a free chiral fermion on a Riemann sphere [1] , or, via bosonisation, as a path 
integral in the theory of a free boson ip with the action: 

S=—l dipdtp+lJ) Jw' 1 + i / Jw+<t Vtf" 1 (J H1 + ...) (2.3) 
4tv J S 2 hJ Ce hJ Ci/e J Cl Y 

with the boson <p normalised so that it takes values in a circle of finite circumference. 
In (2.3) C r denotes a circle \w\ = r, and e [ 0. Similar path integrals were studied 
before in the context of so called conformal matrix models [19], though with very 
different properties of genus expansion [20] . In this paper we shall concentrate on the 
quasiclassical computations in the theory (2.3). 

We start with reminding the basic facts about partitions and free fermions. 

2.1. Partition function and partitions 

The partition X of the size |A| = Ai + A2 + . . . + Xi x is a non-increasing sequence of 
non-negative integers Ai G Z>o (some details and examples are collected in Appendix A): 

A = (Ai > A 2 > A 3 . . . > \t x > Xe x +i = Xt x +2 = • • • = 0) (2.4) 

For our purposes it is convenient to encode the partition A in the so-called profile function 
f\(x), which is a piece- wise linear function, given by: 

f\(x) = \x- a\ + 

00 

+ ^2 \x-a-h(Xi-i + l)\-\x-a-h(Xi-i)\ (2.5) 

i=i 

— I x — a — h(l — i) I + \x — a + hi\ 

The Chern character of the "universal sheaf" ch(£^) (at the fixed point in instanton moduli 
space, characterized by partition A, see [1][2] and references therein) is essentially the 
Fourier transform of the profile function: 

ch(£) = \ J dx fl{x) e ux = e ua ^1 + (1 - e~ uh ) e uh{1 ~ l) (e uhX > - 1) j (2.6) 

hence for the coefficients of its expansion 

00 

MS) = (e* - e"*) ch A (|,^)=E ^ ch fc (a, A) (2.7) 

k=0 
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which enter the formula for the statistical weight of the partition A in the ensemble (1.8), 
one has 

i r 00 

ch fc (a, A) = - / dx f';{x)x k ~ (( a + h ^ " ^ + 1))" " (a + H\ ~ *))*) (2-8) 
Z J »=i 

In the nonabelian case, for the gauge group U(N), the universal sheaf £ splits as a sum of 
N rank one sheaves £/. Accordingly, at the fixed point of the torus action on the moduli 
space of instantons, corresponding to the iV-tuple of partitions A: 

N 00 k 

ch(£) = ( e * -e"*) ^e^ch A(O (0,M = ^ ^ ch fc (a,A) (2.9) 

Z=l fc=0 

with 

JV 

ch fc (a, A) = ^ch fc (az,Az) (2.10) 

For the empty partitions A^ = the Chern character (2.9) reduces to the generating 
function of the vacuum expectation values of the single trace operators tr <E> fc , in the absence 
of quantum corrections. 

Plancherel measure and the profile of the partition 
The Plancherel measure 

j-r Xj-Xj+j-i _ tt (h-iV- TT Xi-Xj+j-i , . 

mX ~U -11 (£ A + Ai -z)! 11 (2 - 11} 

can also be expressed in terms of the profile function f\(x): 

= exp (-± jf ^ cfaicfaa /jfe) /fta*) 7^1 - s 2 )) (2.12) 

where the kernel 7n(x) solves the following difference equation: 

7a (a; + ft) + 7?>0r - ft) - 2 lh {x) = logx 2 (2.13) 
and is given by the asymptotic series: 

7R ( X ) = J_ f F (x) - — F^x) + -FW(x) - — F( 6 )(j) + ...") (2.14) 
mv > h 2 y v ; 12 v ; 24Q v ; 604g v > j v > 



determined by its leading term, or the function: 

F(x) = x 2 (log x - (2.15) 

which enters the perturbative prepotential of pure M = 2 super- Yang-Mills theory. The 
finite difference derivative of 7/1 (x) is related to the V- function: 

1h(x+*;)-1h(x-Vi =log(W(£ + ±Y) (2.16) 



2y " l v 2y ° v 2 , 

In the h — > limit (2.13) just turns into F"(x) = 21og x. The function 7ft(x) is related to 
the generalized Riemann zeta function: 



7»W = -C'(-i) + >gS + ^ 



s=0 



1 I" 00 (it e~ tx 
m'o T* S (e«-l)(e-«-l) (2 ' 17) 



The £'(— 1) and log/i terms in (2.17) are chosen so as to ensure 7r(0) = 0. These corrections 
do not give any contribution to quasiclassical part and can be therefore simply absorbed 
into normalisation (in particular we did not write them in (2.14)). 

2.2. Fermions and Baker- Akhiezer functions 

The gauge theory partition function (1.8), (1.9) can be compactly written as a matrix 
element in the infinite wedge representation of the group GL(oo), in other words, it has a 
free fermion representation [l][3](see also [21]). We shall recall it now in order to motivate 
the introduction of some quasiclassical objects, like the multi- valued functions S(z) and 
$(2), which will be the main tool for solving our problem. As we shall see, they appear 
in the asymptotic expansion of the fermionic one-point functions, or the Baker- Akhiezer 
functions (cf. [22] [23]). 

Free fermions and partitions 

Introduce the free fermion fields, tp(w),ip(w): 



1 



reZ+2 



1 



(2.18) 



r6Z+2 
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where{^ r , ip s } = S rs . The vacuum with the charge M is defined as: 

|M) = |M; 0) = ^_ M+ i V-M+fV'-M+l • • • 
and is annihilated by the fermion harmonics: 

Vv|M) = 0, r > -M, 
Vv|M) = 0, r < -M 

We always use the normal ordering : (...) : with respect to the vacuum |0): 

: ip r ip s : = tfjrtfjs , s < 

: Vv^s : = - ip s ip r , r > 

Let us introduce the W\ +00 algebra operators: 

h k f ~(( / i x fc+lN 



k + 

' 1 



rez+h 



where D = wd w and k > 0. For example, 

Wi = -J ,W 2 = ftL , ... 

i.e. Jo is the zero mode of the U(l) current, 

J(w) = : $(w)ip(w) : = J k w~ k — 

fcez 



J fc = ^ : VvVv+fc : 



while the L Virasoro generator is the zero mode of the stress-energy tensor: 

T(w) = - : $(w)dif>(w) : = ^ L k w~ k (—\ 

kez \ w / 



L k = E r : 



r+fc 



rez+i 
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The importance of the fermions is the relation between the partitions and the excited 
states in the fermion Fock space: 



|M; A) = V-M+I-A 1 V'-M+|-A 2 V'-M+|-A3 • • • V'-M-i+i-Ai • • • ( 2 " 26 ) 



The operators Wfc introduced in (2.22) are diagonal in the basis of "partition" states: 

W fc |M;A> = ^ch fc (ftM,A)|M;A> (2.27) 
with the eigenvalues given by (2.8). In what follows, we shall also use the formula [24] [25]: 

e^jMjO) = ^m A 7r |A| |M; A) (2.28) 



and its direct consequence 



Vw e-^ \m + 1; 0) = h-r-2 y: II • - l M ; A ) ( 2 - 29 ) 

A 71 i=l 1 

where the sum is now taken over all partitions do not containing an eigenvalue, corrse- 
ponding to r-th fermionic mode, which is automatically taken into account vanishing of 
the product factor. The infinite product in (2.29) is actually finite 

fi «-A i + r-i-M _ 1 ^ z-A. + r-1-M 

11 i r( r j- 1 _ m\ 11 i j- t — I _ m ^ ' 



i=i 



Baker- Akhiezer functions 

The sum over partitions (1.8) can be compactly written as a matrix element in the theory 
of free fermions [1] [3] : 

Z(a,t,h) = (M|e-Tre^>° tfcWfc+1 e— \M) = TT^ ^ Lfc >« fc+1 (2.31) 

where a = TiM for M G Z, and the second equality follows from (2.28), (2.27). Let 
r G Z + | and set z = hr. Consider the Baker- Akhiezer function, or the following matrix 
element: 

~ ! Ji ~ i + w J -i 

V(z,a,t,h) = — — (M\e~- ij- r e sL fe >o^ w Hi e — |M + 1) (2.32) 

Zyttj t, fl) 
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Using (2.28), (2.29) it can be expanded into a sum over partitions 

h>0 (2.33) 
• 2^m|(-n _2 ) |A| e^^fe>o fe fe+1 _Q ; — ? 

A i=l Z 

The latter can be obtained from the partition function (2.31) by the shift of times tt — > 

tk ~ S k {h, r): 

*(z,a,t,ft) = ft"-* exp -L £ |^-(( r +i)*+i_ (r . 

" fc>0 

(2.34) 

T(r+i) Z(a,M) 

defined via the generation function 

? E;^= lo «4>tir + s log n^i) (2 ' 35) 

From (2.34) one gets, that asymptotically at Ti — > and far = 2 — > 00 

a, t, ft) ~ exp g( ^ a,t) (2.36) 

with the singularity 

S(z, a ,t) = ^2 f kz k - z{\ogz - 1) + alog^ + . . . (2.37) 

fc>0 

while the nonsingular at z — > 00 terms in (2.37) are expressed, as follows from (2.34), in 
terms of the first derivatives of log Z(a, t, %) w.r.t. t-variables. The essential singularity 
(2.37) contains the Eguchi-Yang term z(\ogz — 1), (cf. [17]), coming from the Gamma 
function logr(^) ~ z(\ogz — 1) in (2.34). The fact that the regular term of the first order 
in z, t\z in (2.37) is accompanied by the z(\ogz — 1) term, of the degree 1 + e, in a sense, 
is the consequence of the asymptotic freedom of the noncommutative U(l) theory (see [3] 
for the discussion of this phenomenon). In the U(N) case we shall see the analogous shift, 
t\ — > t\ — Nlogz, in agreement with the exact coefficient, N, of the beta function of U(N) 
super- Yang-Mills theory [26]. 
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One can also study the fermion matrix element in the coordinate representation: 



%(w t a,t t h) = . 1 „. (M|e~£ e ^fe>o* fcWfe+1 M w ) e^lM + 1) 
Z{a,t,n) 



(2.38) 



Asymptotically, (2.38) gives 



with the quasiclassical phase: 



*n E(w, a, t) 

#(w,a,t,fi) ~ exp ' (2.39) 



E(w,a,t) = S(z,a,t) -z logtu, Q ' ^ = logw (2.40) 

dz 

related to (2.36) by the Legendre transform, while (2.38) can be thought as an integral 
duality transformation [27] for the Baker- Akhiezer functions. Note (for fixed (a, t)): 

= _ Z *E (2.41) 
w 

The nonabelian theory partition function also has a representation in terms of free 
fermions. For the gauge group U(N) one uses N flavours, ip( l \w), I = 1, . . . , N. For the 
special values of the vacuum expectation value of the Higgs scalar ai = h(Mi + pi), 
where Mi e Z are integers, Pi = jj (N + 1 — 2/), the partition function of the nonabelian 
theory can be related to that of the abelian one by means of the procedure of blending of 
partitions. In the language of free fermions this is the following procedure: 

1 N 

V(w-*r) = ^Y^ i - l \w)w- pl (2.42) 

The single free fermion \1/ is produced out of iV free fermions I = 1,...,N living 
on the w-space. The fermion ^ lives on the AT-fold cover of the original space, which in 
quasiclassical limit turns into a spectral curve. 

The spectral curve for a nonabelian theory will be explicitly described below. On a 
small phase space (with only t\ nonvanishing) it coincides with the Seiberg-Witten curve, 
but when higher Casimirs are switched on, it goes beyond the Seiberg-Witten family, 
still being a hyperelliptic curve for the U(N) gauge theory. This curve corresponds to a 
quasiclassical hierarchy [6] of the Toda type, and we shall see that co-ordinate z = hr, 
corresponding to the Baker- Akhiezer functions of the type (2.36) is distinguished from the 
point of view of this hierarchy. 
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2.3. Genus expansion 

The partition function (1.8), (1.9) has the same genus expansion, as in (1.1), but with 
the switched on times t. There are various interpretations of this expansion, which will be 
presented in detail elsewhere [28] . 

a. ) Gauge theory interpretation in four dimensions: the ^-background acts as a "smart 

box" with parameter fi being the inverse volume, prepotential as the extensive part 
of the free energy, the higher terms F g being the finite size effects. 
6.) Topological string interpretation: in the geometrical engineering setup F 5 's are the 
genus g topological string amplitudes on a local Calabi-Yau manifold. 

b. ) Van-der-Vaals gas in one dimensions: the ^-expansion becomes the Mayer diagram 

expansion [28] in this case, 
r.) We can view (2.3) as the quantum field theory and perform the standard Feynman 
diagram expansion around the Gaussian action J g2 dtpdip with the vertices of va- 
lency p weighted with the weight U p ~ 2 , so that the sources § dtpw^ 1 have the weight 
ft -1 , the "mass term" t\ § (dip) 2 having no ^-dependent factor at all, etc. The weight 
of the Feynman diagram 7 with v p {p/) vertices of valency p would be, then 

ftEp^" 2 )^) = ft2(E( 7 )-V( 7 )) = n 2g-2 

where 

£(7) = iZ>p(7) 
p 

is the number of edges, and 

p 

is the number of vertices, and g is the number of loops. In this way one can make 
contact with the string loop expansion. 

2.4. The quasiclassical solution 

The partition functions (1.8) (1.9) should be viewed as the partition functions of an 
instanton gas in the "box" of size TT 2 . We are interested in the free energy of the instanton 
gas, per unit volume: 

F(a, t) = Lim^o h 2 logZ(a, h, t) (2.43) 
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Energy and entropy 

The evaluation of the asymptotics (2.43) is greatly facilitated by the observation that in the 
limit where fi — > with the rest of the parameters kept fixed, the sum over the partitions 
A (or iV-tuples A) is dominated by the single partition of the very large size ~ h~ z [3]. To 
determine this partition and to evaluate (2.43) one can use the geometric representation 
of partitions, Young diagrams. In the limit fi — > the boundary of Young diagram of 
the master partition becomes a continuous curve. This curve, e.g. the profile f\(x), can 
be found by maximizing the weight of the partition A in the sum (1.8). The Plancherel 
measure (2.12) and the Chern characters (2.8) in (1.8) combine together into the following 
energy functional 

E(\, h) = \f dx,dx 2 f»( Xl ) f'i{x 2 ) lh (xi - x 2 ) ~ i / dx %{x)t(x) (2.44) 

4 Jx 1 >x 2 2/1 J 

so that 

Z(a,t,%) =J2 e ~ E{X,H) ( 2 - 45 ) 

A 

Expanding the energy E(X, h) in h (assuming f\(x) — > f(x)): 

h 



where 



£[f] = \f d Xl dx 2 f"(x 1 ) f"(x 2 ) F( Xl - x 2 ) - \ f dx f"(x)t(x) , (2.46) 

we can approximate the partition function as: 

Z(a,t,ft)~exp(-^£[/]) (2.47) 

Indeed, the entropy contribution, i.e. the number of partitions A, whose profiles f\(x) 
differ from f{x) by 0(h), is of the order of 

exp Ql) , 

where L is finite in the % — > limit. Thus in the % — > limit the entropy contribution is 
negligible, compared to that of the energy. 
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The variational problem 

Thus, we arrive at the following conclusion: 



The quasiclassical tau-function or generalised Seiberg-Witten prepotential is the critical 
value of the functional 



F(a,t) = - Ciit f £[f] 



(2.48) 



where we look for the extremum of the energy functional £[/] in the class C 1 (R) of differ- 
entiable functions f(x), such that 



N < f'(x) < N 
f"(z) > 



(2.49) 



For a in appropriate domain of the moduli space M the support of f"{x) is a set of n 
disjoint intervals {li}i=i,..., n along the real axis: 





supp/" = np =1 i h 










x 1 < . 


• < x t < x l+l < x l+l < ■ 





(2.50) 



The moduli a enter the variational problem via the additional constraints: 



1 

ai= 2 


■ / dx xf"(x) , 




k 



(2.51) 



where I;, I = 1, . . . , n is the Z'th connected component of supp/". 

The number n of cuts, and independent moduli ai depend on the gauge group. For 
example, for G = U(N),SU(N), n = N. For G = SO(N), n = 2 [f ] , and for G = 
USp(2N) there are n = 2N + 1 cuts. 

In the G = SU(N) case the Higgs eigenvalues a; obey: J2iLi a i = 0- In the G = 
SO(N) case the Higgs eigenvalues have the form: (ai, a2, . . . , a n /2 5 — a i 5 — a 2, • • • , — a n /2, *), 
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where for odd N, * = and for N even * is absent. In the G = USp(2N) case the Higgs 
eigenvalues have the form: (<zi, 02, ... , ajv, — a±, — . . . , — a at, 0, 0). 

It is easy to notice [29] that our variational problem is very similar to arising in 
the context of multisupport solutions for matrix integrals. However, there are important 
distinctions, caused basically by properties of the functions (2.49), which are extremising 
the functional (2.48). Below we present the solution to the variational problem in the same 
geometric terms, though involving sometimes the multivalued differentials on hyperelliptic 
curves. We shall return to the parallels between our approach and the matrix models, 
when discussing the Eguchi-Yang matrix integral. 
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3. Solution in the U(l) case 

In this chapter we study in detail the extended prepotential in the case of U(l) theory. 
In the Z7 (1) case, the constraint (2.51) can be in standard way taken into account with the 
help of the Lagrange multiplier 

£[f] -> £[f, a]=£[f]-a D [a- 1 -j dxx f"(x)j (3.1) 

The Lagrange multiplier a D naively looks like the zeroth time to, since its contribution 
can be viewed as the shift t(x) — >■ t(x) +a D x. However, in our setup (as often happens for 
similar problems, see e.g. [30] [31]) the variable a is fixed, while a D is varied, thus the dual 
variable a plays the role of the zeroth time. For N > 1 the situation is more complicated 
and will be discussed later on. 

3.1. The limit shape 

The variational equation for the functional (3.1) gives 

-\ J dx f"(x) (x - x) {\og\x - x\ - 1) + t'(x) =a D x E supp/" (3.2) 
It is convenient to introduce the following analytic multivalued function 

S(z) = -\ J dx f"{x) (z - x) {\og{z -x)-l) + t'{z) (3.3) 

on the z-plane with the supp/" removed. It has the following properties: 

1. ) The differential <&(z) = dS/dz is multivalued. However, the differential d<& is already 

well-defined on C, the double cover of the z-plane, which is ramified at the end-points 
of the interval of the support of /". 

2. ) The exponential exp (die $) is therefore single-valued on C. On the cut it is equal to 

one. 

3. ) The equation (3.2) implies that the real part 

UeS(x) = \ (S(x + i0) + S(x - i0)) = a D 

is constant on the cut I = supp/". In order to consider the asymptotic of (3.3) in 
what follows we shall always choose a branch, which is real along the real axis. 
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4.) Asymptotically, as z — > oo, 

OO OO 

fWr - 1\ -I- n 4- V +, ? k - 

kz k dt k 



OO OO r\ sr-" 

S(z) = (log*-l) + alog* + £^ fc -£— — (3-4) 



k=i k=i 

where, according to (2.48), 



dx x k+1 f"(x) , k > (3.5) 



at fc 2(fc+l) 

and the coefficient in front of the z(logz — 1) term is fixed by 

J dx f"{x) = f'(+oo) - f'(-oo) = 2 (3.6) 

5. ) The Legendre transform of S(z), S(w), 

,^ , T , / / / / \ zdz f dx fix) , „. 

dE = -zd§ = d(t' - zt") - — / J —^- (3.7) 

expands near z = oo, w = oo as follows: 

E(z) = -z - - l)t k z k - a \ogz + Y^^rgf (3.8) 

k=2 fe=l ^ fc 

The formulae (3.5) together with the equation 

= f , (3.9) 

which follows from (2.48), identify the generating function (2.48) with the logarithm 
of quasiclassical tau-function, being in the case of the single cut, a tau-function of 
dispersionless Toda chain hierarchy. Note, that the asymptotics (3.4) follows from 
(2.33) (2.34). Note also that the term, linear in z in (3.8), can be viewed as the 
regularised k = 1 term (k — l)tkZ k , where t\ is replaced by the divergent term tq + t±. 
The bare coupling To is a logarithmically divergent function of the energy cutoff [3] . 

6. ) The pair of multi- valued functions (3>(z,a),z) can be viewed as the quasiclassical 

analogues of the pair of Orlov-Shulman and Lax operators [22] . The canonical trans- 
formation, taking the pair (logiu,a) to (<£», z) is generated by the generating function 
S(z,a), or, equivalently, £($, a). 
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3.2. Toda as presentiment 

As we have already discussed, the limit shape is described by the function f(x), 
which differs from \x — a\ on the single interval I = (x~ , x + ). The corresponding integral 
transforms S(z) and suggest to study the double cover C of the z-plane with the 

branch points z = x + and z = x~ . 

The curve C can be equivalently described by the equation 

z = v + A^w + ^j (3.10) 

in C x C*, with x^ = v ± 2A. It is just a copy of a Riemann sphere or CP 1 , with two 
marked points _P±, with z(P±) = oo, w ±1 (P±) = oo. Another form of the curve (3.10) will 
be generalised in the nonabelian case: 

y 2 = (z- v) 2 - 4A 2 = (z - x + )(z - x~) (3.11) 

where y = k(w — to -1 ). The formula (3.3) defines a function with a logarithmic cut and 
asymptotic behaviour (3.4), odd under the involution w <-> — of the curve (3.10). In terms 
of the variable w one can globally write 

S = -S E y(w) + tktik(w) + a logw + a D , (3.12) 

k>0 

with the Eguchi-Yang term 

S E y(w) = + A (w + logw + A(logA - 1) (w - , (3.13) 
the " Hamiltonians" 

VL k {w) = z% - z k _, k>0 (3.14) 
which are the Laurent polynomials in w, and the "zeroth Hamiltonian" 

Q (w) = logw , (3.15) 

all odd under w <-> — . The Eguchi-Yang term (3.13) is fixed by its asymptotics Sey(w) = 
±z(logz — 1) + 0(1) at P± (cf. with (2.37)). Note that one can write the differentials dflt 
in the following suggestive form: 

dVL k = k{z k - 1 y) + - (3.16) 

+ y 
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where the symbol (•••)+ denotes now a polynomial part in z, in the expansion near P + . 
We can think of (3.12) as of the Legendre transform of the Seiberg-Witten differential 
dT, = —z^-. Obviously, on the cut (3.12) gives S\ w=1 = a D . The canonical Toda chain 
times are defined by the residues 

t = resp + dS = —respdS = a (3-17) 

and 

t k = \ res P+ z~ k dS = ~ res P _z~ k dS, k>0 (3.18) 
k k 

From the expansion (3.4) at z — > oo it also immediately follows, that 

|^ = res P+ z k dS = -res P ^z k dS, k > (3.19) 

Notice that (3.4) has no constant in z term, which together with explicit formula (3.12) 
allows to compute a D from a regularized value of the function (3.12) at z — > oo 

« D = s\ w=1 = 

fsr^ k ^ ( 3 - 20 ) 

= lim \2_^t k z - z(\ogz - 1) + alogz + Sey{w) - 2_^tkO, k (w) - alogw I 

2 ^°° \fc>0 A;>0 / 

The consistency condition for (3.19) is ensured by the symmetricity of the second deriva- 
tives 

' = res P+ (z k dQ n ) (3.21) 



dt n dt k 
where we have introduced 



\ n>0 / 

— -Q -d (\o z V 1 

da ° \ da 2 dadt n nz n 

\ n>0 

The functions Q k , k > 0, form a basis in the space of meromorphic functions with poles 
at the points with P±, odd under the w w _1 involution. All time-derivatives here are 
taken at constant z. 
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The expansion (3.22) of the Hamiltonian functions (3.14) expresses the second deriva- 
tives of T in terms of the coefficients of the curve (3.10), e.g. 



, a v A2 + ir 

logz - logA - — 

z z z 



+ 



2A 2 



2vA 2 



z — v — 



+ 



(3.23) 



2 ,2 nk 2s ^K 2 2A 2 (A 2 + 2t> 2 ) 
n 2 = z 2 - (v 2 + 2A 2 ) '- 



+ 



as z — > oo, which gives, in particular, 

da 2 



logA, 



<9 2 .F 
dadti 



and 



a 2 ^ 
at 2 



= 2A 2 = 2 exp 2 



da 2 



(3.24) 



(3.25) 



which becomes the long- wave limit of the To da chain equations after a derivative with 
respect to a is taken: 



2„D 



da 



2^- exp I 2 



da 



D 



da 



(3.26) 



dt\ ~da 

with the Toda co-ordinate a D = The other expansion coefficients (3.23), with the 
help of (3.21), give rise to the Losev-Shatashvili polynomials. In the U(l) case they are 
polynomials in the single variable u = v, whose coefficients depend on A. 

One can now find the dependence of the coefficients of the curve (3.10) on the defor- 
mation parameters t of the microscopic theory by computing dS at the ramification points 
w = ±1, where z = v ± 2A and dz = 0: 



dS 
dlogw 



w=±l 



5> 

fc>0 



dQk 
dlogw 



+ a-vT 2AlogA = 



(3.27) 



w=±l 
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Small phase space 

If tk = for k > 1, solution to (3.27) immediately gives 



v = a , A = e* 1 (3.28) 

and the prepotential 

J= = - aa D + -res P+ (zdS) - \a 2 = - (a 2 h + e 2tl ) (3.29) 
The slope fix) of the corresponding limit shape profile f(x) is given by: 

f\x) = ^arcsin (%^f^j , x e {a - 2e f \ a + 2e l1 ) (3.30) 

where the branch of arcsin (hence the name of the "arcsin law" [32]) is chosen so that 
/' (a ± 2c* 1 ) = ±1. The limit shape itself is given by the function: 

f(x) = - ^A/4e 2tl - (x - a) 2 + (x - a)arcsin ^ *T" )) ' x e (a - 2e f \a + 2e l1 ) 

(3.31) 

and 

f(x) = \x — a|, otherwise. (3.32) 



Turning on £2 

Now let us turn on the Casimir tr $ 3 in the ultraviolet, or in other words consider nonva- 
nishing £i,£2- Then one finds from (3.27) 



v = a _ _L L ^_ m 2 e 2(t 1+ 2t 2 a) j 

logA = *! + 2t 2 a - |L (-m 2 2 e 2 ^ +2t ^) 
where L(t) has an expansion: 



and obeys the functional equation 



(3.33) 



L (t) = ( n £ - t n =t-t 2 + ^t 3 -^ + ... (3.34) 

n=l 



L(£) exp L(£) = t 
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(3.35) 



For the prepotential one gets instead of (3.29) 



T = -aa D + -resp + (zdS) — -^-resp + (z 2 dS) — 



= t(a) + "-e 2t "^ + 2t\ e 4t '» + ^ e 6t '» + 




(3.36) 



where the instanton expansion is governed by the parameter q = e T ° related with the 
bare coupling, which can be introduced by t\ — > t\ + log q. Here tk = for k > 2 and 
t(a) = + t 2 ^-, so that t"(a) = h + 2t 2 a. 

The slope of the limit shape is again given by the trigonometric functions: 



where v and A are the functions (3.33) of ti,t 2 ,a. Hence we see, that the Vershik-Kerov 
"arcsin" law is deformed, for nonvanishing t 2 , by the Wigner semicircle distribution. 

Note that our "rules of the game" are such that the higher order times tk, k > I 
are considered to be the nilpotent parameters. Had we viewed them as the ordinary 
perturbation parameters (and thus faced the possible non-renormalisability concerns), the 
formula (3.37) would have suggested a possibility for the interesting critical behaviour, 
similar to the one found in the two dimensional Yang-Mills theory [33] : indeed, as soon as 
t 2 becomes sufficiently large, \t 2 \ > t 2 , so that 



the slope f'(x) reaches +1 before x reaches the right end of the cut (v — 2A,v + 2A). 
The equation (3.38) is nothing but the convergence radius of the series (3.34) for the 
function L(t). For \t 2 \ > t\ a new cut opens up and the theory goes over to another 
phase, in particular the curve C would grow in genus. Of course, the critical value of t 2 is 
exponentially large ~ A -1 , like the Landau pole. 

3.3. Symplectomorphisms and dispersionless hierarchy 

The evolution of the prepotential jF(a, t) with respect to the times tk is described by 
the Toda chain hierarchy. We start by reminding the Lax formulation of this hierarchy 
[22]. The flows can be described in the Hamiltonian fashion, by interpreting the evolution 




(3.37) 



v-2k<x<v + 2k 



4t* 2 A(a,t 1 ,t* 2 ) = l 



(3.38) 
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to take place on the space of holomorphic functions on V = C* x C. We endow the space 
V with the holomorphic symplectic form: 

vj=— Ada (3.39) 
w 

where w G C*, a G C. The /c'th flow is governed by the Hamiltonian 

tt k (w,a) = z% -z k _ , k>l (3.40) 

where, as in (3.14) (3.10) the ± denote the part, holomorphic at w = or w = oo, respec- 
tively, and z is the Lax operator: 

z(w,a) = v(a) + A(a) (w + , (3.41) 

and 

n = logw (3.42) 
The evolution is encoded in the equations: 
d 



,^z(w,a) = <^z(w,a) , Q k (w,a)j 

dv(a) ( 1 \ dA(a) , A , 

7 +[w+-) = (3.43) 



Thus: 



dt k \ w J dt k 

= A(a) (w - d a n k - (v'(a)w + A' (a) (w 2 + l)) dA 

dfl dtl (3.44) 

= 4 (,(«)A 2 («))' , ^ = A(a) (v(„,f + 2A 2 (a))' 

The first line in (3.44) implies the Toda equation (3.25). Let us discuss the geometry 
behind the equations (3.44) (3.43). They define a set of commuting symplectomorphisms 
of V. Infinitesimally these symplectomorphisms look like: 

Sw = wd a Vtkbtk, = —wd w Vtk8t}~ (3.45) 

Let us denote the finite time symplectomorphism by: 



g (t):V^V, {W(w,a;t), A{w, a; t)) = g{t){w, a) 
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(3.46) 



so that g(0) = id. Now let us take the Lax operator, viewed as a function on V, and let it 
evolve in the geometric way: 

z(w,a;t) = g(t)*z(w, a; 0) = z (W (w, a; t) , A(w,a;t)) = 

v (A (w, a; t)) + A (A (w, a; t)) (w(w, a; t) + — J - —^ j = 

v(a; t) + A(a; t) ( w + — 
\ w 

the last equality being a consequence of the asymptotic behaviour of the Hamiltonian 
vector fields (3.45) as w — > and w — > oo. 

We can also change the canonical variables, from (it?, a) to (z, $), such that: 

dw 

— A da = dz A (3.48) 

10 

By substituting (3.41) into (3.48) we can express $ via it; and a (of course, (3.48) determines 
$ up to an addition of a function of z): 

$(w,a)= [ da z{w,a) = v(a)+A(a) (w + — \ (3.49) 

Finally, the symplectomorphism g(t) can be described with the help of a generating func- 
tion (see [34] for the definitions): 

dS = $dz + logw da (3.50) 

where we view S as a function of (z,a). We can think of (-2,$) as of the functions of 
W, A, i.e. they coordinatise the image of g(t). Thus S is a function of z, a, t. We claim it 
coincides with (3.3). 
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4. The nonabelian case 



4.1. The energy and surface tension 

In the nonabelian case, G = U(N),N > 1, similarly to (3.1), one has to add to the 
functional S[f] several Lagrange multipliers, or the so-called surface tension term: 

=*[/]-£ a? U-k [ dxxf"(x)) = 

(4-1) 

= £[f] + \a D J xf"{x)dx + J dx<7(/')-J>aF 



where 

N 



" JV 



and a is a piece- wise linear function, which vanishes for t outside the interval (—N,N), 
and equals ±(af - a D )(t + N - 21 + 2) + Emiifam - a D ), for t G (-AT + 2Z - 2, -AT + 2/), 
for I = 1, . . . , AT, so that: 



V | xr(x)dx + | efe <t(/') = | £ of | x /" (*) rfx (4.3) 
Define, as before: 

S(z) = -\J dx f"(x)(z-x)(\og(z-x)-l)+t'(z) (4.4) 

and 

= = - \ J dx f"(x)\og(z -x) + t"(z) (4.5) 

4.2. The curve 

The variational equation, following from (4.1), can be written as: 

S( x ) = \ (S(x + iO) + S(x - iO)) = af , for x G I; (4.6) 

which implies: 

$(x + iO) + $(x - iO) = , for x G 1/ (4.7) 
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Together with the obvious property 

S(x ± iO) = af =F — / dx f"(x) (x - x) , for xeli (4.8) 
it implies that the differential 

= S"{z)dz = ~ Q dx dz + t"'dz (4.9) 

is well-defined on the double cover C of the z-plane, which is ramified at the points x^: 

N 

y 2 = H(z-x+)(z-x^) (4.10) 
i=i 

Let us denote by P± the two preimages of the point z = oo. The differential <i$ is odd 
under the involution y i— > — y. As z — > oo, on one of the sheets of C, say near P + , the 
differential <i$ has the following asymptotics: 

fc>i 2 Z fc>i 

These two properties allow to write: 

d$ = Y^ s k z k — = ^k{k- l)t fc dO fc _i - NdQ - 2%i ^ dwj (4.12) 

fc>0 ^ fc>l j = l 

where the differentials dflk are fixed by their asymptotics at z — > oo 

r ±/cz fc - 1 ^ + 0(z" 2 ), £;>0 
dO fc ~ <^ (4.13) 
I ±f + 0(z- 2 ), k = 

and vanishing A-periods 

dO fc = 0, fc>0, Vi = l,...,iV-l 

Ai 

while doji, i = 1, . . . , N — 1 are canonical holomorphic differentials normalized to the bi- 
cycles, surrounding the first N — 1 cuts. 

The coefficients Sk and the ramifications points xf 1 are to be determined from the 
conditions: 

— / z~ k+1 d$ = ±k(k - l)t k , A; = 2,3,... 
' " Jp± (4.14) 

res p ± d§ = 
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2ni 



directly following from (4.11), 

-^nf A M= ll f " {x)dX = * ( //( ^ +) " //( ^ _) ) = 1 > 



1 / (4 - 15) 
_ J) rf$ = 0, 1 = 1,.. .,N 



B, 



If s(z) = J2k>o s k zk is a polynomial of power N + K — 2 (in the case of nonvanishing times 
ti, . . . , tx up to the K-th order), its higher K coefficients are fixed by leading asymptotic 
and the residue at infinity via (4.14), while the rest N — 1 coefficients can be determined 
from fixing the A-periods by the upper line in (4.15). Now it leaves only 2N branch points 
{xj} of the curve (4.10) to be fixed from the rest of (4.15) (vanishing of the S-periods), 
and the Seiberg-Witten periods of the "dual differential" 

<f cE = --L / Z d$ = ai , l = l,...,N (4.16) 
2m J A[ 2m J A[ 

Together with normalization of $>(xjj) = at the largest positive ramification point, 
suggested by integral representation (4.5), the second conditions of (4.15) can be equiva- 
lently written as 2N real conditions 

&e®{xf) = 0, j = l,...,2N (4.17) 

The condition (4.17) ensure that under variation at constant z 

5 (dS) = 5 (<&dz) ~ holomorphic 

modulo the exact terms. 
Since under rescaling: 



Z i— » /J,Z, |/H /j, N y 

xf h-> jixf, ai h-> [iai, t k i-> ^~ k tk 



(4.18) 



the function 3>(z) acquires logarithmic corrections, the t\ dependence can be fixed by 

£i = res P+ (z'^dz) (4.19) 
with the appropriately chosen branch of $ at the point P+. 
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The Lagrange multipliers 

af = §^, 1 = 1 JV 

oai 

can be computed by a standard trick. From (4. 6), (4. 8) it is clear that 



af - af 

or 



dT _ 1 
dai 2 j B . 
and still 

1 N 1 N XT 

D _ \ ^ D _ \ ^ 

a ~ N^ ttl ~ N^da t 

i=l i=l 



J J dS = ^<J) dS (4.20) 
dS, z = l,...,iV-l (4.21) 



For the other time-derivatives one can write 

dT 

dt k 



res P+ (z fc dS) = - r ^—res P+ (,2 fc+1 d$) (4.22) 



4.3. Small phase space and Seiberg-Witten curves 

Now suppose: 

2 

X 

The conditions (4.11) (4.15) imply that w = exp (— <fr) is a meromorphic function on C with 
the only pole at P+, where it behaves as: 

wocz N + ... , (4.23) 

(recall that z(P+) = oo) and which transforms under the y h- > — y involution aswn m; -1 , 
(for the switched on higher times, such function would acquire an essential singularity at 
the points where z = oo, and the constraints (4.15) cannot be resolved algebraically). 
These conditions fix w uniquely, up to a few parameters: 

/ 1 \ N 

lw + - = ]J(z - Vl ) = P N {z) (4.24) 
^ w ' i=i 

i.e. 

dZ . AT f 1 



d$ = -P' N (z)-, y = A N ( 
y \ 



w 

w 



The coefficients of Pn(z) are determined by the periods (4.15): 

1 



zd$ = ai, l = l,...,N (4.25) 



2ni 

Finally, the parameter A is related to the time t\ via: 

A N = e fl (4.26) 
In this way one easily recovers the results of [3] . 
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4.4. Quasiclassical flows 

Let us now come back to the general formula (1.7) and show how it can be derived 
from our quasiclassical solution on large phase space, or for nonvanishing microscopic times. 
Taking an extra derivative of (4.22) one gets formula (3.21), which can be re- written on 
generic non-abelian curve (4.10) as 

' " V ~ = Tes P+ (z m dQ n ) = res P+0P+ (z(P) n z(P') m W(P, P')) (4.27) 



dt n dt 



where we have introduced the bi-differential W(P, P') = dpdp>\ogE(P, P'), E(P, P') being 
the prime form [35], with the only second order pole at diagonal and vanishing A periods. 
In the inverse co-ordinates z = z(P) and z' = z(P') near the point P + with z(P + ) = oo 
it has expansion 

1 ' k>0 

The bi-differential W{P, P') can be related with the Szego kernel [35] 

S e (P,P')S- e (P, P') = W(P,P') + du} i (P)du} j (P , )^-log e (O|r) (4.28) 

which, for a half-integer characteristics e = — e, has an explicit expression on hyperelliptic 
curve (4.10) 

S(zz >)- Ue(z) + Ue(z')VdzU 

2y/U e (z)U e (z>) z-z> 



with 

JV 

u e (z) = n 

j 



z - x p + 
J_ 

Z ~ X e~ 



(4.30) 



Here {x e ±} is a partition of the ramification points of (4.10) into two sets, corresponding 

3 

to a given characteristic e. For example, on a small phase space, when (4.10) turns into 
the Seiberg-Witten curve (4.24), there is a distinguished partition e = E with 
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Substitution of (4.28) into (4.27) gives 

d 

- res P+ (z n du z ) res P+ (z m du 3 ) - — log tf e (0|r) = (4.32) 

OTij 

f \ d 2 T d 2 F d 
C nrn [x e ±J Q a .Qf n ddjdtm drij ^ ^ e ^ T ^ 

where for the Losev-Shatashvili polynomials one gets from (4.29) 

r ( \ 1 ( fl ge(g) U°WY \ A A A (A ^ 

C nm (x ef ) = -resp +8P+ [1 + + — j dzdz j (4.33) 

If calculated on the small phase space, where all = 0, with /c > 1, and the particular 
choice of the characteristic (4.31), which correspond to the deformation (1.7) around the 
Seiberg-Witten prepotential, residues in (4.33) vanish for n,m < N, and one gets 

d 2 T du n+1 du m+1 d 

log ve{0\t), n,m<N (4.34) 



dt n dt m da,i da,j drij 

with 



u ™ = = - res ^+ z 



l -H v ?= l -^® n ) ( 4 - 35 ) 



dt n -i n + V y J n jr[ n 

justifying (1.7). 

Our derivation of the renormalisation group equation almost repeats, or even simpli- 
fies, the derivation from [8]. The main difference is that the formalism of quasiclassical 
hierarchy is developed now at fixed hyperelliptic co-ordinate z, while in [8] the role of 
such distinguished co-ordinate was played by the co-ordinate on base torus w. In the last 
case the quasiclassical flows were not deforming the geometry of the Seiberg-Witten curves 
(only the generating differential had to be replaced), but the Losev-Shatashvili polynomi- 
als (the analogs of (4.33)) were never vanishing, being instead related to the generalised 
Kontsevich model, or topological Landau-Ginzburg models [36]. Now we see, that alge- 
braic w has to be replaced by the transcendental <E>. The proper choice of local co-ordinate 
z is suggested by the microscopic theory, where it is already encoded in .Fuv = tr t(<E>), 
see (1.6), and therefore it appears to be distinguished in the effective functional (2.48). 
We can conclude therefore, that the quasiclassical hierarchies arise as adequate language 
for the effective low-energy theories, but to specify the details (the choice of proper basis, 
local co-ordinates etc) one needs to turn directly to a microscopic theory. 
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4.5. Instanton expansion with higher times 

Let us now test our results in the non-abelian case against explicit calculations of the 
first instanton corrections to the prepotential. We shall see immediately that extracting 
even the first term in the instanton expansion of the Krichever tau-function is a non- 
trivial task, and the direct evaluation using the sum over partitions here is much more 
effective. The true application of the quasiclassical tau-function is the analysis of the non- 
perturbative effects such as the emergence of the massless BPS particles at some points 
of the moduli space of vacua, the surfaces of marginal stability and so on. However, 
when we study the extended Seiberg-Witten theory, the "new" directions coming from the 
"times" tk which couple to the higher Casimirs, are to be viewed as nilpotent (much like 
the generalised moduli space of the topological string of B type on Calabi-Yau manifold 
[37]). We should not attempt to give a physical meaning to the singularities which might 
occur at finite values of t^s. Note that adding matter to the M = 2 theory effectively 
switches on the times tk- For example, if we add Nf fundamental multiplets, we get 
tk ~ fc^r m J 2k - I* 1 the theory with "real" fundamental multiplets the prepotential has 
singularities when the charged matter becomes massless, i.e. near a ~ m/, which is near 
infinity for small t. 

Let us start with the perturbative limit with the higher times switched on. It is 
characterized by the degenerate differential (4.12) 



where the position of poles of degenerate holomorphic differentials coincide with the per- 
turbative values of the Seiberg-Witten periods 




(4.36) 



a t = -Tes Vl (zd® ) = v t , l = l,...,N 



(4.37) 



Integrating (4.36), 



N 




(4.38) 



i=i 



one gets the perturbative generating differential dSo = <&odz, satisfying 



ddS 



dz 



1 = 1, ...,N 



dai 



z — 
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ddSo k _-i 
— — = kz k L dz, k > 
ot k 

and, in the rational case, it is easy to write the function 

N 



S (x) = t'(x) - ^(x - a t ) (log(x - ai) - 1) 



1=1 

Equations 

define the perturbative prepotential with switched on higher Casimirs to be 

JV 



= S ( aj ) (4.39) 



^o = ^t(a ; ) + ^^F(a ; -a m ) (4.40) 



2 

1=1 m^l 



to be the sum of the ultraviolet prepotential JF UV = tr t(<3>) and the perturbative M = 2 
prepotential, with the " Seiberg- Witten" function defined in (2.15). 

Consider now vanishing S-period (4.15) of the perturbative differential (4.36) 



Jx+ 



rf$ = 



where x^ = cij ± 2 yj qSj + 0(q 2 ) are positions of the branching points of the curve (4.10) 
in the vicinity of perturbative rational curve. Equation & (x^) = $>o(x~) m ^ ne nrs ^ 
nonvanishing order in parameter of instanton expansion q gives 

e 2t"(a;) 

5, = = ^, 1 = 1,. ..,N (4.41) 

where the numeric coefficient is fixed, say, from comparison with the Seiberg- Witten curve 
(4.24) on a small phase space. 

Take now for simplicity the U(2) gauge group. Then for the equation of the curve 
(4.10) up to the second order in q one has 

y 2 = II {(z - at) 2 - 4qSt + . . .) = 

i=i,2 (4.42) 
= (z- ai ) 2 (z - a 2 ) 2 - 4q (S 1 (z - a 2 ) 2 + S 2 (z - ai ) 2 ) + 0(q 2 ) 

and 

dz dz ( ( S\ S 2 



y (z - a 1 )(z - a 2 ) \ \{z - en) 2 (z - a 2 ) 2 ■ 
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One can also expand therefore (4.12) as cM> = <i$ + qd§\ + 0(q ) with 



.//// s. 2z - ai - a 2 \ f Sx S 2 



dg 1 = 2 t'"(z)- " - - _^_ + 



(z - ai){z - a 2 ) J \(z- ai ) 2 (z-a 2 ) 2 / 
2(5i - 5 2 ) 1 ' J 



(ai - a 2 )(z - ai)(z - a 2 ) 



normalized to res ai <i$i = res a2 <i$i = (we restrict ourselves now for t"'(z) = 2t 2 with 
only nonvanishing t\,t 2 ^ 0). Computation of residues of (4.43) at infinity gives 



- -^resoo (z 2 d^ 1 ) = ai8ai + a 2 5a 2 + Si + S 2 

- Trresoo (z 3 d$i) = a\8ai + a 2 8a 2 + 2 (aiSi + a 2 S 2 ) 

Taking into account that (4.22) gives at linear order in q 



(4.44) 



vesoo {z k+1 d^ x ) = ^- 



= £>f"W§£ (4.45) 



where the first term in the r.h.s. arises as a contribution from the perturbative prepotential 
(4.40) due to instantonic renormalization of the relation between the coefficients of the 
curve (4.10) and the Seiberg-Witten periods as ai — > ai + qSai with 

8a i = (2*2 J #i = res ai {zd^ x ) 

) ai \ a \ (4.46) 

Sa 2 = \ 2t 2 H ] ^2 = res a2 (zd&t) 

V ai - a 2 / 

Comparing (4.45) with (4.44) and taking into account (4.46), one gets T\ = \ Si, what 
coincides with direct instantonic calculus. Higher coefficients of the instantonic expansion 
of the quasiclassical tau-function can be obtained in a similar way. 
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5. Eguchi-Yang picture 

In 1994, T. Eguchi and S. Yang have proposed [17] the following matrix model: 



(t) = I VM exp iVtr ( -2M (logM - 1) + ^ t k M k J (5.1) 
J V k>0 J 



for the description of the type A CP 1 topological string, in the stationary sector of Gromov- 
Witten theory (some generalisations for other target spaces were studied in [38]). Their 
proposal is more general, it includes the full phase space of the sigma model, i.e. the full 
set of gravitational descendants of the cohomology of CP 1 . 

One can easily relate the matrix model (5.1) to our variational problem (3.2). We shall 
show that the proposal [17] is correct as far as the genus zero part of the topological string is 
concerned, if one views (5.1) as the integral over the supermatrices of (N + \N-) x (N + \N-) 
size. 

Indeed, let us study the saddle point approximation to the matrix integral (5.1), in 
the large N limit. As usual, we assume the eigen-values to form a distribution, with the 
continuous density of the eigenvalues: 

p( x ) = L {trS (x-M)} (5.2) 

The action of the matrix model (5.1) together with the contribution of the measure (the 
Vandermonde determinant) combine together nicely to the following effective action func- 
tional of p(x): 



V ef f[p] = 2 dx 1 dx 2 p(x 1 )p(x 2 )log(x 1 -x 2 )+ / dxp(x) -2x(log(x) - 1) + J^t 

Jx 1 >x 2 J V b^n 



k>0 

(5.3) 

In the case of bosonic matrix model the density function p(x) is a non-negative function 
on the real line which has a compact support. The obvious constraint for (5.2) 



dxp(x) = 1 (5.4) 
can be enforced with the help of a Lagrange multiplier: 



VefT[/9]^VefT[p]-a D (^J dx P (x) - 1 
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The saddle point equation reads now: 

2 / dxp(x)\og\x — x\ — 2x(log(x) — 1) + y^t fc x fc = a D , x G suppp (5.5) 

^ fc>0 

Note that once we shift x — > x — a, make an identification: 

/'(x) = sgn(x) - (5.6) 

and rescale the times tk and a D by an irrelevant here factor of 2*, the equation (5.5) 
becomes equivalent to (3.2). It can be shown by integrating (3.2) by parts and choosing 
the real branch of the logarithm function. However, the "density" p(x) which corresponds 
to the saddle point of (3.1) is not non-negative. For example, on the small phase space, 
for a = 0, 

p(x) = sgn(x) - -arcsin (7^7) (5.7) 

which has a jump at x = 0, where it goes from being positive to being negative. Note also 
that for (5.7) obviously: 

J dx p(x) = (5.8) 

(which is also consistent with (5.6) and (2.51)). 

What could this negative density of eigenvalues possibly mean? The hint comes from 
(5.8). Since the integral of the density of eigenvalues is proportional to the trace of unit 
matrix, one can conclude that the "matrix model" corresponding to the original model of 
random partitions (1.8) must have vanishing trace of the unit matrix. The natural venue 
for such matrices is the supermatrix model (see, e.g. [39]), where the traces are replaced 
by supertraces. For the (N + \N-) x (N + \N-) matrices the supertrace of a unit operator is 
equal to N + — iV_, which can be both positive and negative. 

M= (c d) ' strM = trA - trD (5.9) 

where A and D are iV + x iV_|_ and iV_ x iV_ size hermitian matrices with bosonic entries, 
respectively; B and C are fermion-valued iV+ x iV_ and iV_ x iV_|_ dimensional complex 



£ The factor of 2 in (5.1) in front of the Eguchi-Yang term is a consequence of the fact that 
matrix integral is rather related to the correlation function bilinear in fermions, cf. e.g. with 
[19]. In particular, this is related to the fact that integration measure contains the square of the 
Vandermonde determinant. 



38 



matrices, B = C" 1 ". The only subtlety in the definition of the measure DM is the division 
by the volume of the gauge group, which is done by the proper treatment of ghosts (both 
bosonic and fermionic). The invariance of the measure and the action (which is the same 
as (5.1) with the replacement Ntr — > j^str, with the new "supermatrix Planck constant" 
h) is such that both A and D can be diagonalised: 

A->diag(a^,...,:r^ + ), D -> diag(xf , . . . , x^_) (5.10) 

We can view (5.10) as a gauge-fixing condition. The corresponding Fadeev-Popov deter- 
minant equals (note again its similarity with the (square of the) free fermion correlator): 

Yll<i<j<N + ( X i ~ X j ) lll<i<j<7V_ ( x i ~ x j ) 



lll<i<Ar + Yll<j<N^ ( x t X j ) 

Now one can write for the density of eigenvalues (5.10): 



(5.11) 



P (x) = h[Y t s(x- xt) - £ - x i) ( 5 - 12 ) 

\i=l 1=1 J 

and consider it in the large iV limit, i.e. when h — > as N± — > oo, so that hN± are finite. 
The logarithm of jacobian (5.11) plus the Eguchi-Yang action are still given by p-V e ff [p], 
defined in (5.3), but the supermatrix density (5.12) is automatically subjected to (5.8) 
upon iV + = iV_, which can be re-expressed in terms of times of the Toda hierarchy. For 
example, on the small phase space: 

hN± -> 2^^e tl , (5.13) 

TV 

while in general: 

hN± -> ^(t) -a\ (5.14) 

It would be nice to investigate our supermatrix model further, perhaps extending the sad- 
dle point correspondence to the exact relation between the h and (or /i)-expansions. 
Perhaps it is worthwhile to study more general supermatrix models, with iV+ ^ N-. It is 
also interesting to find out whether there is a relation of our supermatrix model and the 
Dijkgraaf-Vafa matrix models, computing the superpotentials of M = 1 supersymmetric 
gauge theories, as suggested in [1][3]. In this paper we have pointed out new features of 
the matrix models possibly relevant for the M = 2 physics: the non-positivity of eigen- 
value densities, which could be related to the block-matrix decompositions arising in the 
perturbative matrix model calculations (cf. e.g. with [40], [41] where such decomposition 
was performed for bosonic one- and two- matrix models). 
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6. Conclusions and future directions 



In this paper we have solved the long-standing problem of extending the prepotential 
of the effective M = 2 theory to the " large phase space" , obtained by turning on the chiral 
perturbations 



where t is a gauge-invariant function. We have seen that for the single-trace perturbations 
the resulting prepotential coincides with the Krichever tau-function of the quasiclassical 
hierarchy. For the U(l) theory this is the dispersionless Toda chain hierarchy. For the 
U(N) theory this more general quasiclassical or universal Whitham hierarchy (in the sense 
of [6]) of the Toda-chain type (i.e. it corresponds to hyperelliptic curve with two marked 
points P±). 

The derivation of the paper was based on a particular identification of the operators 
(6.1) with some cohomology classes of the moduli spaces Aik of gauge instantons, and 
their continuation on the partial compactifications Aik- For the gauge theory with gauge 
group U(N) the compactification Aik is the moduli space of rank N torsion free sheaves 
on X = CP 2 with the fixed trivialisation at the fixed complex line D = CP^ (this is also 
a moduli space of noncommutative instantons [9]). The operator (6.1) of the form: 



is identified with the k-th component of the Chern character ch(£) of the universal sheaf 



One can imagine other definitions of the ultraviolet observables (6.1). This discussion is 
parallel to the discussion of the 2-observables in the two dimensional sigma models [10]. 
The integrable structure which we find with our prescription is very suggestive. Of course, 
the arguments [42] which show that the moduli space of vacua of Af = 2 supersymmetric 
gauge theory is a base of an algebraic integrable system, based on electric-magnetic duality, 
do not apply per se to the large phase space, spanned by the couplings However, one 
can think, at least formally, of them as coming from the degeneration of the standard 
moduli space of vacua of the theory with larger gauge group, more matter and so on. 

The appearance of quasiclassical hierarchies in the Seiberg-Witten theory was antici- 
pated long time ago [4] . The program of studying the " large phase space" of the Seiberg- 
Witten theory was put forward in [7] and some proposals about corresponding integrable 




(6.1) 




(6.2) 
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hierarchies were made in [4] [8] [43] . The subtlety of all these proposals is the contact term 
ambiguity. In a sense it is a standard problem of reconstructing the ultraviolet theory 
given the infrared one. We have shown in this paper, that the quasiclassical hierarchy, 
in the basis of time-variables proposed by microscopic instanton theory, solves the ambi- 
guity problem exactly in the way anticipated in [7]. For small k,l the Losev-Shatashvili 
polynomials can be also studied using the blowup techniques [7] [44] . 

One should keep in mind the richness of physics which is related to the random par- 
tition models of which we studied the thermodynamic limits: our results can be easily ex- 
tended to the five dimensional supersymmetric gauge theories (which means going beyond 
the small phase space analysis of [23]) where they become the statements about the disper- 
sionless hierarchies governing the melting of the crystals with appropriate symmetries (thus 
generalizing the work of [45]); the same fermion correlators and their dispersionless limits 
show up in the physics of one dimensional electrons [46] ("physics of quantum wires"). 
Finally, one may study more seriously the analogues of Douglas-Kazakov phase transition 
on the large phase space, and also incorporate the multi-trace ultraviolet perturbations. 
We hope to return to this discussion elsewhere. 
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Appendix A. Partitions and Chern characters 

Partitions and Young diagrams (the standard reference is [24]). 

For each partition A (2.4), its size \X\ = Ai + A2 + . . . + Xg x is a number of boxes in the 
so-called Young diagram Y(X), which is a geometric way to represent partition. The length 
£\ is the number of columns. Traditionally one draws Young diagrams as the collection of 
rows of boxes, of the lengths Ai, A2, etc. Here we draw an example of the Young diagram, 
for the partition A = (5, 2) with £\ = 2, |A| = 7: 



and its profile function f\(x) defined in (2.5), 



Profile of the partition 5,2 




Fig.l 



For partition A one introduces the dual partition A', whose Young diagram Y(X') is 
obtained by exchanging the rows and columns of Y(X). For example, the partition, 
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dual to (5, 2), 



, is the partition 
(2,2,1,1,1): 



It is convenient to introduce a coordinate system, on the Young diagram: 

1 < j < A, < l x , , 



1< % < A' < £ x . 



(A.l) 



The coordinate % labels rows, from top to bottom, while the coordinate j labels columns, 
from left to right. For example, the partition A = (5, 2) gets coordinatised in the following 

way: 



(1,1) 


(1,2) 


(1,3) 


(1,4) 


(1,5) 


(2,1) 


(2,2) 





For the box with the coordinates the hook-length is defined as: = Xi + X'j — i 

j + 1. In our example A = (5, 2): 



tyl.l) = 6 


h(i,2) = 5 


h(i, 3 ) = 3 


^(1,4) = 2 


^(1,5) = 1 


h( 2 ,i) = 2 


h(2,2) = 1 





Finally, the Plancherel measure, which appears in (1.8), (2.11) is just the product of the 
hook-length's 



1 



—nf =n 



{h-iV- -q A, - Xj + j - % _ -q A, - Xj +j - i ^ A2 ^ 



J ~ * 



i<j 



J - * 



□ e > y^A + A.-z)!^.^ 
In our example: 

1 l!-0! 4 4 78 9 10 11 12 345 



m (5 , 2 ) = 



1-1-2-2-3-5-6 6! -2! 1 12345 6 7 123 
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Chern characters. 

The Chern character of the partition A is a generating function for the symmetric functions 
of the eigen-values \ — i: 

00 1 „Mu 00 1 

ch A (M,«)= Ve u(M+ r^) = „ e ^ + y e <M +1 -i) {e uX i _ 1) = 

f— f e 2 — e 2 ^— f 

i=i i=i 

1 / M 2 1 \ 

= - + M + M (--- + |A|j+ (A.3) 

+ y(^-^ + 2Af|A| + £A,(A,-2 t + l)) + ... 

the second equality being valid when 5teu > 0. In the gauge theory the main object is 
the Chern character of the "universal sheaf" £, which is related to (A.3) in the following 
simple way: 



, . / hu hu \ / Qj \ 

ch(S) = [e 2 - e 2 J c h A huj 



e ua ( 1 + (1 - e -^)^ e ^(i-*)( e ^ _ ^ 
i=i 

OO 



(A.4) 



fc=0 

It is the components chfc(a, A) of Chern character (A.4) that enter the formula for the 
statistical weight of the partition A in the ensemble (1.8). For example: 

ch (a, A) = 1 
chi(a, A) = a 

ch 2 (a,A) = a 2 + 2ft 2 |A| (A - 5) 
ch 3 (a, A) = a 3 + 6ft 2 a|A| + 3ft 3 X t (X t + 1 - 2i) 

i 

One can also write for the Chern characters 

00 

ch fc (a, \) = a k + ((a + H\ ~ » + !)) fc - (a + ft(A* - i)) k - 
i=i 

- (a + h(l - i)) k + (a - M) k ) = (A.6) 

OO 

= ((« + ^ - i + l)) k - (a + ft(A, - z)) fc ) 
i=i 

which is useful to make link with the fermionic formalism. 
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